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Let /at and £Z represent, respectively, the free energy per spin and the inverse 
spin-spin correlation length of the critical Ising model on a N x oo lattice, with 
/jV — > foo as N — > oo. We obtain analytic expressions for Ofc and 6^ in the expansions: 
N(f N ~ foo) = EfeLi afc/iV 2 ^ 1 and ^ = E£=i b k /N 2k ~ l for square, honeycomb, 
and plane-triangular lattices, and find that bk/au = (2 2k — l)/(2 2fc_1 — 1) for all 
of these lattices, i.e. the amplitude ratio b^/a^ is universal. We also obtain similar 
result for a critical quantum spin chain and find that such result could be understood 
from a perturbated conformal field theory. 
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Experimental data, analytic and simulational studies of phase transition models, and 
renormalization group (RG) theory suggest that critical systems can be grouped into uni- 
versality classes so that the systems in the same class have the same set of critical exponents 
|];|2|||. RG theory has also been used to propose that critical systems of the same universal- 
ity class could have universal finite-size scaling functions (UFSSF's) and universal amplitude 
ratios 0||,|| and some analytic and numerical calculations of critical systems have supported 
the idea of universal amplitude ratios @,§|Q]- Using Monte Carlo methods H and choosing 
appropriate aspect ratios for lattices of critical systems, Hu, et al. have found UFSSF's 
for percolation and Ising models 0. Despite the success of RG theory and Monte Carlo 
simulations, it is valuable to have more analytic results which could widen or deepen our 
understanding of the universality of critical systems. In this Letter we present exact calcu- 
lations for a set of universal amplitude ratios for the Ising model on square (sq), honeycomb 
(he), and plane-triangular (pt) lattices [0,|],|§ and for a quantum spin chain ITOj , which is in 
the universality class of two-dimensional (2D) Ising models. As far as we know, no previous 
RG arguments, analytic calculations, or numerical studies predict the existence of this whole 
set of universal amplitude ratios. 

Let /jv and £^ x represent, respectively, the free energy per spin and the inverse spin-spin 
correlation length of the Ising model [0,^,0 on a N x oo lattice with periodic boundary 
conditions, with f N — > f^ as N — > oo. In this letter, we obtain analytic equations for a k 
and bk in the expansions: 



oo t 

£n =22 jy2fc-l' ( 2 ) 

fc=l 

for sq, he, and pt lattices, and find that 

b k /a k = (2 2k - l)/(2 2fe " 1 - 1), (3) 

for all of these lattices, i.e. the amplitude ratio a k /bk is universal. We also obtain similar 
expansions for the critical ground state energy E Q and the critical first energy gap (Ei—E ) of 
a quantum spin chain [nj, which are, respectively, the quantum analogies of the free energy 



and inverse spin-spin correlation length for the Ising model, and find that the amplitude 
ratios have the same values. We could physically understand such result from a perturbated 
conformal field theory. 

Consider an Ising ferromagnet on an iV x M lattice with periodic boundary conditions 
(i.e. a torus). The Hamiltonian of the system is 

(3H = -JY, s iSj , (4) 

where (3 = (fc^T) -1 , the Ising spins Sj = ±1 are located at the sites of the lattice and the 
summation goes over all nearest-neighbor pairs of the lattice. We consider a transfer matrix 



acting along the M direction |H|,[12|,[13| . If A and Ai are the largest and the second-largest 



eigenvalues of the transfer matrix, in the limit M — ► oo the free energy per spin, fjy, and 
the inverse longitudinal spin-spin correlation length, £^ , are 

f N =-^\nA and ^ = ^\n(A /A 1 ). (5) 

Here £ is a geometric factor which is 1, 2/v^3 and l/v^3 for sq, he, and pt lattices, respec- 
tively || . Exact expressions for eigenvalues Ao and Ai are available for all lattices under 
consideration: sq |T|, [n| , |i2| , [T3lJ , he [@|J, and pt ||. 



We start from the Ising model on the sq lattice. Onsager has obtained expres- 
sions for all eigenvalues of the transfer matrix. The two leading eigenvalues are A = 
(2 sinh 2J) N ' 2 exp(|Ef=o 1 72r + i), Ai = (2sinh2J) 7V / 2 exp (i Ef=i 72r), where 7fc is implic- 
itly given by cosh7fc = cosh2Jcoth2J — cos (kir/N). At the critical point J c of the sq lattice 
Ising model, where J c = | In (1 + y/2), one then obtains 7& = 2tjj sq (|^). Here 



ip sq {x) = In ( sin 2 + y 1 + sin x\. (6) 



Then the critical free energy /# and critical spin-spin correlation length £jy of Eq. (|5|) can 
be written as 

1 1 N - x 

1 JV-1 

^ = 77 Z!(72r+l-72r). (8) 

/ r=0 

It is readily seen from Eqs. (^|), (0) and (||) that ^ x and Nf^ have odd parity as a function 
of A" -1 . Therefore, in the following expansions of £^ and A^/at as a function of AT -1 , we 
keep only odd terms. 

To write /# and ^ in the form of Eqs. ([[]) and (@), we must evaluate Eqs. (0) and 
(|) asymptotically. These sums can be handled by using the Euler-Maclaurin summation 



formula EM. After a straightforward calculation, we have 
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Here _B 2 fc are the Bernoulli numbers and il)f q k ^ = (d 2h 1 ip sq (x)/dx 2k 1 ) _ ; b 2 = 7r 3 /96 has 
been computed previously by Derrida and Seze [O]. 

For the Ising model on the honeycomb (he) lattice, Husimi and Syozi || found that 
A = (2sinh2J) Jv exp (71 + 73 + ... +7at-i) and 



Ai = (2sinh2J) Ar exp (^j + 72 + ••• + 7v-2 + \in)i where the 7 r are given by cosli7 r = 
cosh 2 J cosh 2 J* - sin 2 f - cos f (sinh 2 2Jsinh 2 2J* - sin 2 f f' 2 . Here J* is defined by 
(cosh 2 J — 1) (cosh 2 J* — 1) = 1, so that one has J* = J at the critical point (J c = 
\ In (2 + a/3)) and one then obtains 7 r = 2-0^ (f§0, where 



1pkc(x) 



\nU(x) + ^A 2 (x)-l\, (11) 



with A(x) = (a/8 + cos 2 2x — cos 2a;) /2. Using the Euler-Maclaurin summation formula, 
we can write the free energy f^ and the inverse spin-spin correlation length £^7 for the he 
lattice as 



2k-l 



Mff f \ V^ v / 3J? 2fc (2 2fc - 1 - 1) / it y- ' (■*_!) 
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+ TT7T Uu? +-' ( 12 ) 



t _i _ ^ v oil 2H z ~ x ; / n \ /(2fc-i) 

4iV " fe tl (2jfe)! UivJ ;V ' 



7T 3 / 71 \ 5 93 / 7T \ 9 

~ 4iV ~ 10 \3n) + 10\3n) + '"' ^ ' 

where iphc = \d 2k ~ l iphc{x) I 'dx 2 ^ 1 ) _ . 

For the case of the pt lattice, we note that one can use the star-triangle transformation 
to transform the he to the pt lattice [|llj] . The amplitudes of the A^ -3 and A^ ~ 7 correction 



terms are identically zero for the he and the pt lattices. 

Above results for the sq, he, and pt lattices can be summarized as 

where a = 1, ( = 1, tjj(x) = 4> S q(x) for the sq lattice, a = 1, ( = 2/\/3, ip(x) = iphc(x) 
for the he lattice, and a = 2, ( = l/v3, ip( x ) — 4>tr(x) = 4>hc(x) for the pt lattice. The 
ratios of the nonvanishing amplitudes of the JV~( 2fc-1 ) correction terms in the spin-spin 
correlation length and the free energy expansion, i.e. bk/dk, axe the same for all three 
lattices under consideration |TJ|]. Thus we have established Eq. @. 

To check whether Eq. (|3|) is still valid for other models in the Ising universality class, 
we proceed to study a quantum spin model on a one- dimensional lattice of A" sites with 
periodic boundary conditions, whose Hamiltonian is [|l(J 

\ N 1 N 

H = -- E < ~ 77 E [(1 + 7K + i< + (1 - 7K + i<] , (16) 

z 7 n =l 4 1 n=l 



where a x , a y and a z are the Pauli matrices. The phase diagram is well known [|T^]. For all 
7 (0 < 7 < 1), there is a critical point at A c = 1, which falls into the two-dimensional Ising 
universality class. The inverse correlation length ^~ 1 is given by the difference in eigenvalues 
Ei — E Q of the Hamiltonian H. In particular, the first energy gap gives the inverse spin-spin 
correlation length ^f 1 and the second energy gap is the inverse energy-energy correlation 
length £^~ 1 . By expanding the exact solution of Eq. (0), Henkel |17]] has obtained several 
finite-size correction terms to the ground-state energy Eq, to the first (E\ — Eq) and second 
(E 2 — Eq) energy gaps. We have extended the calculations to arbitrary order and found that 
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E 2 — Eq = 22 
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where , # 2fc x) = (gP 1 ip q (x)/dx 



„2k-l 



x=0 



ip q (x) 



sin 



non-universal number a = - Jq ip q (x)dx = 21 + arccos 7/(7a/1 — 7 2 ) 
of amplitudes for (E x — Eq) and (— E ) also satisfy Eq. ([3]). Equations 



1 x — (1 — I/7 2 ) sin x, and cvq is an 

jix. Thus, the ratios 
17| and (|IJ) implies 



also that the ratios r& of amplitudes for (E 2 — Eq) 
by 



4A- 



Tk 



(2 2k ~ 



l)B 



and (—Eq) are 7-independent and given 



(20) 



2k 



It is of interest to compare this finding with other results. The exact and numerical 
estimates |Tj| of the subdominant correction amplitudes for the sq, he and pt lattices are 
presented in Table I, which shows that the numerical values obtained by Queiroz 



arc 



very close to our exact results. On the basis of conformal invariance, the asymptotic finite- 
size scaling behavior of the critical free energy and the inverse correlation length is found to 



be 19 



N 



lim N\f N -f 



en 



(21) 



lim JVff 1 = lim N(E t - E ) = 2nxi, 



N^oo 



N 



(22) 



where c is the conformal anomaly number and x« is the scaling dimension of the ith scaling 
field. For the 2D Ising model, we have c = 1/2, x\ = rj/2 = 1/8 and x% — 1 and the leading 



terms of Eqs. flT4]) and (|15|) for all of the sq, he, and pt lattices are consistent with Eqs. 
(PI) and (|^), i.e. ai and &i are universal. Equations ( ETJJ) and (|22"D implies immediately 
that their ratio is also universal, namely 

lim T777 ^T = r l and lim T777 ^T = f l' ( 23 ) 

N^oo N(f N - /oo) *-«, N(f N - /oo) 

where ri = 12xi/c and fi = \2x-ijc. For the 2D Ising universality class we have ri = 3, f\ = 
24, which is consistent with Eqs. ([3]) and ( p0| ) for the case k = 1. 

The corrections to Eqs. ( [2~lf ) and (|^) can be calculated by the means of a perturbated 
conformal field theory f20| , [2~If . In general, any lattice Hamiltonian will contain correction 
terms to the critical Hamiltonian H r 



r N/2 

H = H e + '£g p J <f> p (v)dv, (24) 



where g p is a non- universal constant and (f> p (v) is a perturbative conformal field. Below 
we will consider the case with only one perturbative conformal field, say 4>i(v). Then the 
eigenvalues of H are 

rN/2 

E n = E n>c + g t < n\<h(v)\n >dv + . . . , (25) 

J-N/2 

where E nc are the critical eigenvalues of H . The matrix element < n\4>i{v)\n > can be com- 
puted in terms of the universal structure constants (C n i n ) of the operator product expansion 
|2T)|| : < n\<pi(v)\n >= (2tt/N) Xi Cu n , where x\ is the scaling dimension of the conformal field 
(f>i(v). The correlation lengths (£~ x = E n — E Q ) and the ground-state energy (E ) can be 
written as 

C 1 = ^xn + 27rgi(C nln - C m ) (jp) ' + . . . , (26) 

E Q = E Q)C + 2n gi Com (jf) +■■■■ (27) 



Equations ( |2q ) and (|27| ) show that while the amplitude of correction to scaling terms are 
not universal, ratios of them are. For the 2D Ising model, one finds |22| that the leading 
finite-size corrections (1/iV 3 ) can be described by the Hamiltonian given by Eq. (|24] ) with a 
single perturbative conformal field <pi(v) = L 2 _ 2 {v) + L 2 _ 2 {v) with scaling dimension x\ = 4 
. The universal structure constants C212, Cm and Cqiq can be obtained from the matrix 
element < n\L 2 _ 2 (v) + L 2 _ 2 {v)\n >, which have already been computed by Reinicke |23| 



(C212 = 1729/5760, C m = -7/720 and C m = 49/5760). Equations © and (|7]) implies 
that the ratios of first-order corrections amplitudes for (E n — E ) and (— E ) is universal 
and equal to (C i — C n i n )/C i , which is consistent with Eq. (^) and Eq. ( |20"D for the cases 
n = i,A? = 2 [(Coio - C m )/C m = 15/7] and n = 2, k = 2 [(C olo - C 2l2 )/C m = -240/7] 
respectively. Compare the amplitudes of the iV -3 correction terms for the Ising model and 
the quantum spin chain with the general results of Eqs. fl2HD and (E7|) one can find that 
9i — (3/7 2 — 4)/567r for the quantum spin chain and gi = — 1/287T for the Ising model on 
the sq lattice. For the Ising model on the he and pt lattices we find that gi = 0, which 



indicate that at least two perturbative conformal fields are necessary to generate all finite- 
size correction terms. Further work has to be done to possibly evaluate exactly all finite-size 
correction terms from perturbative conformal field theory. 

The results of this Letter inspire several problems for further studies: (i) On the basis of 
perturbated conformal field theory, can one find other universal amplitude ratios? (ii) How 
do such amplitudes behave in other models, for example in the three-state Potts model? 
(iii) For the critical Ising model on a large iVxMsq lattice (M/N is a finite number), we 
have obtained expansions in N~ l for the free energy, the internal energy, and the specific 
heat [53]. It is of interest to extend such expansions to inverse spin-spin correlation length 



and to he and pt lattices, and to study whether the amplitude ratios are also universal. 
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TABLES 

TABLE I. Comparison of exact (Eqs. (^), (|To|) , ( |i~2| ) - (|l5[)) and numerical |lj|] values for 
subdominant finite-size corrections terms in free energy and inverse spin-spin correlation length 
expansion. 
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Table I. 



Izmailian and Hu 





square 


honeycomb 


triangular 


exact 


numerical 


exact 


numerical 


exact 


numerical 


a 2 


0.15072... 


0.150730(2) 





< 10" 6 





< 10" 6 


h 


0.322982... 


0.322987(6) 





<10" 8 





<10" 8 


a 3 


0.39213... 


0.385(1) 


-0.18590... 


-0.1865(10) 


-0.01161... 


-0.01165(5) 


h 


0.79692... 


0.790(1) 


-0.37780... 


-0.3777(2) 


-0.02361... 


-0.02360(1) 


0,4 


3.13146... 














K 


6.28759... 














a 5 


48.9925... 




7.03535... 




0.02748... 




h 


98.0809... 




14.0844... 




0.05501... 





